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We iiave measured tlie current-voltage characteristics of small-capacitance single Josephson junc- 
tions at low temperatures (T < 0.04 K), where the strength of the coupling between the single junc- 
tion and the electromagnetic environment was controlled with one- dimensional arrays of dc SQUlDs. 
We have clearly observed Coulomb blockade of Cooper-pair tunneling and even a region of nega- 
tive differential resistance, when the zero-bias resistance of the SQUID arrays is much higher than 
the quantum resistance h/e^ ~ 26 kQ. The negative differential resistance is evidence of coherent 
single-Cooper-pair tunneling in the single Josephson junction. 



PACS numbers: 73.23.Hk, 73.40. Gk, 74.50.-hr 

Small-capacitance superconducting tunnel junctions 
provide a novel system for studying the interplay be- 
tween the Josephson phase and the charge on the junc- 
tion electrode. These are quantum mechanically conju- 
gate variables, and their behavior is influenced by dissi- 
pation, such as discrete tunneling of quasi-particles, cou- 
pling to an electromagnetic environment, etc. The sim- 
plest and the most fundamental example is the single 
junction. Current-voltage {I-V) characteristics of single 
junctions have been the subject of extensive theoretical 
investigations Experimentally, however, the obser- 

vation of charging effects such as "Coulomb blockade" 
has been considered to be extremely difficult in single 
junctions because a high-impedance environment is nec- 
essary, and special care should be taken with the mea- 
surement leads . For this reason, thin- film resistors ||] 
and tunnel-j unction arrays were employed for the 
leads, and an increase of differential resistance around 
V = was reported. 

We use one-dimensional arrays of dc superconducting 
quantum interference devices (SQUIDs) for the leads. 
The advantage of this SQUID configuration is that the 
effective impedance of array can be varied in situ by ap- 
plying an external magnetic field perpendicular to the 
substrate. Thus, the zero-bias resistance of the SQUID 
arrays at low temperatures can be varied over several 
orders of magnitude. This phenomenon has been exten- 
sively studied in terms of the superconductor-insulator 
transition The single junction in our samples, on 
the other hand, does not have a SQUID configuration, 
and therefore its parameters are practically independent 
of the external magnetic field. This enables us to study 
the same single junction in different environments. This 
type of experiment has not yet been reported to the best 
of our knowledge. We show that the I-V curve of the 



single junction is indeed sensitive to the state of the en- 
vironment. Furthermore, we can induce a transition to a 
Coulomb blockade of the single junction when the zero- 
bias resistance of the SQUID arrays is much higher than 
the quantum resistance Rk = h/e^ « 26 kSl. 

In addition to Coulomb blockade, we have clearly 
observed a region of negative differential resistance in 
the I-V curve. Negative differential resistance has 
been reported in one-dimensional arrays |^ and in two- 
dimensional arrays however, clear observation in sin- 
gle junctions has not yet been reported in the literature. 
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FIG. 1. A scanning electron micrograph of a sample. 



In the context of the theory of current-biased single 
Josephson junctions [0j2j, the negative differential resis- 
tance appears as a result of coherent tunneling of sin- 
gle Cooper pairs. The onset of the coherent tunneling 
is characterized by the local voltage maximum in the 
low-current part of the I-V curve, or blockade voltage 
Vfc. The value of Vf, is suppressed as the ratio Ej/Ec 
is increased. Here, Ej is the Josephson energy and 
Ec = e^/2C, where C is the capacitance of the single 
junction. We show that the Ej/Eq dependence of our 
measured Vh is consistent with the theoretical prediction. 
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A scanning electron micrograph of an AI/AI2O3/AI 
sample is shown in Fig. |l|. A single junction with the 
area of 0.1x0.1 is in the center of Fig. ^ On each 
side of the single junction there are two leads enabling 
four-point measurements of the single junction. A part of 
each lead close to the single junction consists of an array 
of dc SQUIDs. The area of each junction in the leads is 
0.3x0.1 fim^ and the effective area of the SQUID loop is 
0.7x0.2 fim^. All of the samples have the same configu- 
ration, except for the number N of junction pairs in each 
lead. In Fig. |l|, = 17 is shown, however, we measured 
the samples with larger N, up to 225. The tunnel junc- 
tions are fabricated on a Si02 substrate using electron- 
beam lithography and a double-angle-evaporation tech- 
nique 0. The oxidation conditions, or the thickness of 
the AI2O3 layer, determines the tunnel resistance of the 
junctions. The samples are characterized by the normal- 
state resistance i?„ of the single junction, the normal- 
state resistance r„ per junction pair of the SQUID arrays, 
and N. 
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FIG. 2. Zero-bias resistance of two SQUID-array leads 
connected in series for a sample with r„ — 1.1 kfi and 
N = 225. 

The samples were measured in a ^He-^He dilution re- 
frigerator mainly at the base temperature (0.02 — 0.04 K). 
The temperature was determined by measuring the resis- 
tance of a ruthenium-oxide thermometer fixed at the mix- 
ing chamber. To measure the resistance of the thermome- 
ter, we used an ac resistance bridge (RV-Elekroniikka, 
AVS-47). The samples were placed inside a copper RF- 
tight box which was thermally connected to the mixing 
chamber. All the leads entering the RF-tight box were 
low-pass filtered by 1 m of Thermocoax cable We 
measured the I-V curve of the single junction in a four- 
point configuration, where the potential difference was 
measured through one pair of SQUID-array leads with 
a high-input-impedance instrumentation amplifier based 
on two operational amplifiers (Burr-Brown OPAlllBM) 
and a preamplifier [Stanford Research Systems (SRS) 
SR560]. Through the other pair of SQUID-array leads, 
the bias was applied, and the current was measured with 
a current preamplifier (SRS SR570). When the voltage 



drop at the SQUID arrays was much larger than that 
at the single junction, the single junction was practically 
current biased. The SQUID arrays could be measured 
in a two-point configuration (same current and voltage 
leads) on the same side of the single junction. Note that 
the two arrays are connected in series and that current 
does not fiow through the single junction. In order to 
obtain i?„ and r„, we measured at T = 2 — 4 K (above 
the superconducting transition temperature of Al). 

The effective Josephson energy between adjacent 
islands in the SQUID arrays is proportional to 
I cos(7ri?^/$o)l> where B is the external magnetic field 
applied perpendicular to the substrate, A is the effective 
area of the SQUID loop, and $0 = h/2e = 2 x lO-^^ Wb 
is the superconducting flux quantum. Figure ^ shows 
the zero-bias resistance ro as a function of _B, of two 
SQUID arrays on the same side of the single junction. 
The arrays have r„ — 1.1 kfi and N = 225. In our 
samples with A = 0.14 /im^, BA/^q becomes 1/2 at 
B = 7.4 mT. Thus, tq should have the first local max- 
imum at 7.4 mT, as observed in Fig. || The value of 
tq at the local maximum is more than 10^ times larger 
than that at B = 0. This means that we can tune the 
electromagnetic environment for the single junction over 
a wide range. Henceforth we indicate the magnetic field 
as the frustration / = BA/^q. 
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FIG. 3. Current-voltage characteristics of the single junc- 
tion for the sample with R„ = 17 kf2, r„ — 1.4 kfl, and 
— 65. From top to bottom, the frustration / = BA/^o is 
0, 0.22, 0.41, 0.45, and 0.49, respectively. The origin of the 
current axis is offset for each curve for clarity. The vertical 
broken lines and the boxes with dotted lines represent the su- 
perconducting energy gap ±2A and the region that will be 
closed up in Fig. W, respectively. 



In Figs. |3|-|5|, we show some results on a sample with 
R„ = 17 kfl, r„ = 1.4 kn, and N = 65. The I-V curves 
for the single junction at several frustrations are shown 
in Fig. 0. The vertical broken lines represent the su- 
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perconducting energy gap ±2A, which is ±0.34 meV for 
Al The features within the gap are dependent on 

/ especially in the vicinity oi V = 0. A close-up of the 
region denoted by the boxes with dotted lines is shown 
in Fig. ^. As f is varied, the single-junction I-V curve 
shown in Fig. Q develops a Coulomb blockade. We reem- 
phasize that the Josephson energy of the single junction 
is independent of /, because it does not have a SQUID 
configuration and the field f^o/A applied here is much 
smaller than the critical field for Al films (~ 0.1 T). The 
electromagnetic environment for the single junction (the 
SQUID array), however, is strongly varied with /. Fig- 
ures ^ and 1^ demonstrate that the low-bias region of the 
single-junction I-V curve is indeed sensitive to the envi- 
ronment. 
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FIG. 4. Current-voltage (I-V) curves of the single junc- 
tion for the same sample as in Fig. ^ From top left to bottom 
right, the frustration / is increased from 0.43 to 0.49 in steps 
of 0.01. The origin of each curve is offset for clarity. For 
the labeled curves, the I-V characteristics of the leads at the 
same / are shown in Fig. ^ 

The I-V curves of the two SQUID-array leads con- 
nected in series at / = 0.43 (ro 0.61 Mfl), 0.46 
(ro 3.2 MQ), and 0.49 (ro 43 Mn) are shown in 
Fig. |. The I-V curves of the leads are nonlinear, and in 
general the SQUID array cannot be described by a liner 
impedance model However, we may characterize the 
environment by their vq. Coulomb blockade is visible 
only when tq 3> Rk, which is consistent with the the- 
oretical conditions for the clear observation of Coulomb 
blockade in single junctions pT[ |. For an arbitrary linear 
environment characterized by Ze{uj), 'Re[Zf,{u>)] ^ Rk is 



required for the Coulomb blockade of single-electron tun- 
neling and Re[Ze(tj)] ^ Rk/^ for that of Cooper-pair 
tunneling [|ll[. However, as seen in Fig. ||, our SQUID 
arrays are a nonlinear environment. It is interesting to 
note that at / = 0.46 (labeled "b"), the I-V curve of the 
leads is still "Josephson-like" (differential resistance is 
lower around y = 0), while that of the single junction is 
already "Coulomb-blockade-like". This feature becomes 
more distinct in samples with larger N . 

The region of negative differential resistance seen in 
Fig. ^ when Coulomb blockade is well developed, is re- 
lated to coherent tunneling of single Cooper pairs accord- 
ing to the theory j^j^ of a current-biased single Joseph- 
son junction in an environment with sufficiently high 
impedance. For sufficiently low current, the electrical 
conduction of single Josephson junctions is governed by 
stochastic quasi-particle tunneling, and the I-V charac- 
teristic is highly resistive. As the current is increased, 
coherent single-Cooper-pair tunneling, or "Bloch oscil- 
lation" dominates, decreasing the mean voltage. As a 
result, the I-V curve has a region of negative differential 
resistance, or "back-bending" in the low-current part. 
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FIG. 5. Current-voltage curves of the two SQUID-array 
leads connected in series for the same sample as in Figs H 
and ^. From a to c, the frustration is 0.43, 0.46, and 0.49, 
respectively. 

Theoretically, the value of the local voltage maximum, 
or blockade voltage Vb, is a function of the ratio EjjEc-, 
and given by 



Vbr 

as T ^ 0, where 



0.25 e/C ioYEjlEc<.\, 
6o/e for Ej/Ec > 1, 



v3/4 



Ec) 



1/2 



(1) 



(2) 



is the half width of the lowest energy band [|. 
When quasi-particle tunneling is neglected, Vb becomes 
larger ||l|], 
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e/C for E.,/Ec < 1, 
7r(5o/e for Ej/Ec > 1. 



(3) 



The calculation of Vb for arbitrary Ej jEc has also been 
done for the case of no quasi-particle tunneling jl^ . 

In Fig. 1^, we compare the measured Vf, with the above 
predictions. We used Ej = /iA/8e^i?„, and for Ec we 
estimated C from the junction area. A value of specific 
capacitance = 45 ± 5 fF/ /im^ |14| , which was obtained 
for the junctions with 3 x 28 /inr and 7 x 54 /im^, has 
been frequently employed Uncertainty in Cs, how- 

ever, seems to be much larger when the junction area is 
on the order of 0.01 /zm^ or smaller. 

In the following experimental literature on small tun- 
nel junctions, C was estimated from the offset voltage in 
the normal-state 1-V curve. Fulton and Dolan measured 
samples with three junctions that share a common elec- 
trode, and obtained 0.20-0.23 fF for (0.03 ±0.01 /im)^ x 
3 III], i.e., Cs = 42-192 fF//zm2. Geerligs et al. reported 
Cs « 110 fF I [ivc? for two-dimensional (190 x 60) junction 
arrays with the areas of 0.01 or 0.04 [ivc? More re- 
cently, Penttila et al. studied resistively shunted single 
Josephson junctions with the area of 0.15 x 0.15 /im^ [|6| . 
The estimated C of their eight samples ranged between 
0.8 and 6.6 fF, or 36 - 293 fF/^xm^. 
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FIG. 6. Blockade voltage Vb divided by e/C as a func- 
tion of Ej/Ec- From top to bottom, the curves represent 
the theoretical predictions of Refs. [^, [^, and respec- 
tively. For the calculation of the solid symbols (lower data set) 
45 fF //im^ is assumed, and for the open symbols (upper data 
set) 200 £F //zm^. The solid boxes and the open boxes repre- 
sent the samples with the nominal junction area of 0.01 ^rr? 
in Ref. |l. 

Motivated by the above survey, in Fig. ^ we plot the 
experimental Vb vs Ej/Eq, where two different values 
of Cs have been used to determine C . The lower data 
set uses 45 fF//im^ and the upper data set 200 fF j yuxr? . 
The data from Ref. |^ for the samples with the same 
nominal junction area (0.1 x 0.1 /zm^) as in this work are 
also plotted. The measured Vb agrees with the theoreti- 
cal predictions qualitatively. We also find a quantitative 
agreement when Cs ~ 200 fF j [ivc? is taken for 0.01 /im^. 



This apparently large junction capacitance may be partly 
explained by distributed capacitance of the SQUID ar- 
rays. 

In summary, we have studied the current-voltage char- 
acteristics of single Josephson junctions biased with 
SQUID arrays. We have demonstrated that the Coulomb 
blockade in the single junction is indeed sensitive to the 
state of the environment. The region of negative dif- 
ferential resistance has been clearly observed in the 1-V 
curve of the single junction, when the zero-bias resistance 
of the SQUID arrays is much higher than Rk-, which is 
clear evidence of coherent single- Cooper-pair tunneling 
in the single Josephson junction. 
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